Universality of QCD traveling-waves with running coupling 
beyond leading logarithmic accuracy* 

GUILLAUME BEUF''' 

Institut de Physique Thcoriquc, 
CEA, IPhT, F-91191 Gif-sur-Yvette, France 
CNRS, URA 2306 



We discuss the solutions of QCD evolution equations with saturation 
in the high energy limit. We present a general argument showing that, 
in the running coupling case, the Next-to-Leading-Logarithmic (NLL) and 
higher order terms are irrelevant for the universal asymptotic features of 
the solutions. 



1. Introduction 

The BFKL equation pQ is known to provide an incomplete description of 
the hadronic or nuclear collisions in the the high energy limit of QCD. That 
evolution equation relies indeed on the assumption of a dilute partonic con- 
tent of incoming particles, but it inevitably leads to denser and denser ones. 
It also violates the Froissart bound and therefore unitarity. At high density, 
and hence at very high energy, coherent collective effects modify the BFKL 
equation by reducing the emission of additional soft gluons. That mecha- 
nism of gluon saturation tames the BFKL growth of the cross sections, and 
is expected to restore unitarity. Several evolution equations with saturation 
have been derived, like the B-JIMWLK [2 [3] or the BK [21 0] equations, 
implementing nonlinear effects on top of the BFKL equation. However, the 
result discussed here is independent of the precise saturation mechanism at 
work. 

The features of gluon saturation are well known at the leading logarith- 
mic accuracy (LL) with fixed QCD coupling. The saturation scale Q S (Y), 
which is the typical momentum scale for the onset of nonlinear effects, is 
related to the color correlation length in the transverse plane. Hence, that 
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scale plays the role of an infrared cut-off for the gluons radiation. By anal- 
ogy with the Fisher and Kolmogorov-Petrovsky-Piskounov equation, it has 
been understood [5] that gluon saturation also back-reacts onto the dilute 
linear regime. By linearity, the solutions of the fixed coupling BFKL equa- 
tion can be decomposed as sum of wave solutions, interpreting formally the 
rapidity Y as time and L = log/c^/A 2 as space position for the wave (kx 
being the gluon's transverse momentum). If a saturation mechanism is at 
work, one of these waves solutions, the critical one, is selected dynamically 
in a L interval which is growing with Y. Therefore, the solution loses mem- 
ory of its initial condition during the evolution, even in the linear domain. 
The asymptotic shape of the solution and the asymptotic evolution of the 
saturation scale are then universal. That selection mechanism provides a 
natural explanation for the geometric scaling [6] observed in the HERA 
data, as the critical wave solution has precisely that scaling property. 

One should also consider the effect of running QCD coupling and of 
higher logarithmic orders, in order to do reliable studies. A first step in that 
direction consists in replacing by hand the fixed coupling in the LL equation 
by a running coupling, taken for example at a scale equal to the parent gluon 
transverse momentum kx- In that case, the mechanism of selection of the 
critical wave solution seems to hold. However, the running coupling leads 
to a quite different universal asymptotic behavior of the solution [5j [7] (see 
[8] for a review). The evolution is slower, and in particular the saturation 
scale behaves as \og Q 2 S {Y) oc VY, instead of Y in the fixed coupling case. 

On the other hand, one can introduce NLL contributions to the kernel, 
but keeping the coupling fixed, so that the linear part of the evolution 
equation writes formally 

dyN(L, Y) = a [xlU-Ol) + a XNL L (-d L ) + 0(a 2 )] N(L, Y) , (1) 

where a = N c a s /ir, and N is the Fourier conjugate of the dipole target 
amplitude. In that case, the asymptotic behavior of the solution is similar 
to the one in the fixed coupling LL case, except that the value of the critical 
parameters are changed, and depends on the NLL kernel eigenvalues [9]. 



2. Saturation with running coupling and a NLL kernel 

Let us now discuss the most complete case, with both running coupling 
and higher orders. The result presented here was first found in Ref.|10j 
and Ref . [TT] . It is possible to resum some of the NLL and higher order 
contributions by letting the coupling run, for example at the scale kr, i-e. 
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taking a(k^) = 1/bL. Then, Eq.(pQ) rewrites 



dyN(L,Y) 



bL 



XLL{ - dL) + ^\- dL) + 0{(^) 2 



N(L,Y). 



(2) 

Using the w-expansion method of Ref.[12j, one replaces Eq.Q by the effec- 
tive equation 



d Y N(L,Y) = -L Xeff (-d L ,dy)N(L,Y). 



with a kernel whose eigenvalues are 



Xeff(j,u) = XLL{i) 1 + 



Xll(7) 



+ 



(3) 



(4) 



In the brackets in the previous expression, the terms of re-th order to n are 
supposed to have poles of order n at most, at 7 = and 7 = 1. Hence, that 
w-expansion is safe only when u <C 7, 1 — 7- Assuming that condition to 
be fulfilled, if one tries to solve Eq.([3]) in the large Y limit using Laplace 
transform, one gets the position of the saddle-point uj s ~ Y~ x l 2 . Thus, for 
large enough Y (in practice several rapidity units), the w-expansion will be 
safe, and we will have logiV oc \/Y and logQ^(y) oc y/Y '. Expanding the 
kernel of Eq. ([3]) around u> = and changing variable Y into VY, one gets 

^5^iV(L,y) = Xeff(-d L ,0) N(L,Y) 

^fXeff{-d L ,0)d V y + O(d^ 



N(L,Y) ,(5) 



where Xeff stands for the derivative with respect to u. Thanks to Eq.Q, the 
first line in Eq.([5]) is the LL BFKL equation with running coupling, and the 
second line contains the higher order terms. The relevant regime determin- 
ing the universal asymptotic properties of the solution corresponds to L oc 
\[Y 3> 1, and more precisely to r ~ F 1 / 6 > 1, with r = \og{kl/Q 2 a {Y)) be- 
ing the geometric scaling variable. Solving Eq.([5]) perturbatively at large Y, 
one finds that the three leading orders, corresponding to 0(Y°), 0(Y^ 1 ^ 6 ) 
and C(y _1//3 ), features only terms from the first line of Eq.([5]). The lead- 
ing order give the dispersion relation between 7 and uj. In the case of the 
critical wave solution, which will be the relevant asymptotic solution in the 
presence of gluon saturation, the second order selects 7 = j c , solution of 
Xll{1c) = 7c Xi,i(7c)- It determines the leading term of log-/V(L,y) and 
of log Qg (Y). Their next-to-leading terms are given by the 0(Y~ 1 / 3 ) order 
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terms in Eq.([5]). In N(L, Y), it corresponds to the Airy diffusion factor. Up 
to that order, the NLL kernel has not been relevant. 

Therefore, we have shown that in the running coupling case, the uni- 
versal asymptotic solution in the linear regime derived in Refs.[5j[7] is valid 
not only for LL evolution equations with saturation but also for NLL ones 
(or even all orders ones). NLL terms in the kernel can modify the solution 
only at subleading orders not calculated in Refs.0 [5]. For example, the 
third term in the large Y expansion of log Q^(Y), which will be of order 
C^y- 1 / 6 ), will be NLL dependant. 

Saturation and running coupling are thus essential ingredients to un- 
derstand the high energy limit of the partonic content of hadrons, whereas 
other NLL or higher order contribution are not. One should however notice 
that the behavior of the solutions of saturation equations is more subtle in 
the running coupling case, and not fully understood [13], by contrast to the 
fixed coupling case. 
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